Superexponential droplet fractalization as a hierarchical formation of dissipative 

compactons 
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We study the dynamics of a thin film over a substrate heated from below in a framework of a 
strongly nonlinear one-dimensional Cahn-Hillard equation. The evolution leads to a fractalization 
into smaller and smaller scales. We demonstrate that a primitive element in the appearing hierar- 
chical structure is a dissipative compacton. Both direct simulations and the analysis of a self-similar 
solution show that the compactons appear at superexponentially decreasing scales, what means 
vanishing dimension of the fractal. 
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Introduction. — A vast number of intriguing pattern- 
formation phenomena can be described with high-order 
nonlinear diffusion equations of Cahn-Hillard type. Since 
their introduction [1] , these equations have been success- 
fully applied to a great variety of natural and techno- 
logical processes such as phase separation in binary mix- 
tures, alloys, glasses, and polymer solutions (see, e.g., 
surveys [2]), topology transitions in a Hele-Shaw cell [3], 
dynamics of layered systems [4] , thin films [5] , competi- 
tion and exclusion of biological groups [6] , and aggrega- 
tion of aphids on leaves [7]. In the thin film context, nu- 
merical studies of an amplitude equation of Cahn-Hillard 
type [8, 9] have evidenced film rupture leading to the for- 
mation of a cascade of "drops" and "fractal-like finger- 
ing" [10] comprising the gaps or "dry spots" [9] between 
the drops. These findings have been supported by direct 
simulations of the Navier-Stokes equations [12]. 

The goal of this paper is to describe this cascade as a 
hierarchical formation of dissipative compactons. Com- 
pacton is a well-known compact (i.e. with finite support) 
traveling-wave solution, which emerges in conservative 
systems with nonlinear dispersion [11]. Its stationary 
analogue with compact support appears in dissipative 
systems with nonlinear dissipation and, therefore, can 
be referred to as a stationary "Dissipative Compacton" 
(DC). Below we demonstrate that a DC is a primitive 
element mediating the formation of hierarchical fractal 
structure, and characterize the fractal properties of this 
structure quantitatively. 

We focus on a one-dimensional Cahn-Hilliard equa- 
tion describing dissipative evolution of a conserved field 
h{x, t) 



h t + [f(h)h x + g{h)h x 







(1) 



In the context of the dynamics of a thin film over a 
substrate heated from below, this equation describes a 
surface-tension-driven convection (see, e.g., Eq. (4) in 



Ref. [8]), where h is the local thickness and 



/ = -Boh 3 
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g = h 3 



(2) 



Dimensionless Bond (Bo), Biot (B), and Marangoni (M) 
numbers determine the levels of the gravity, of the heat 
flux at the free surface, and of the convective flow, re- 
spectively. Although function / here has a rather com- 
plex form, for h — > one can set / ps 0.5BM h 2 . This 
approximation holds also for moderate values of h, pro- 
vided the gravity can be neglected, Bo —> 0, the heat 
transfer at the free surface is poor (B small) while the 
thcrmocapillary effect is strong (M large). The function 
g(h), which can be referred to as "mobility", is conven- 
tionally non-negative, g(h) > 0, as this prevents against 
the fast growth of the short-wave perturbations. 

Assuming the limiting form of /, after an appropriate 
rescaling of the time, the field, and the coordinate we 
arrive at our basic equation 

h t + (h 2 h x + h 3 h xxx ) x =0. (3) 

Noteworthy, Eq. (3) is invariant under the scaling 

h^p 2 h, x -> px, t -> p~ 2 t , (4) 

meaning that a thinner film evolves slower. 

Steady state and its stability. — We start our analysis 
by considering positive stationary solutions h = H(x) 
of Eq. (3). Looking for symmetric patterns, after one 
integration we obtain 



HH'" + H' = 



(5) 



with primes denoting d/dx. Equation (5) admits a com- 
pact solution H(x) in the form of a DC or a "touchdown 
steady state" [13], nonvanishing for |x| < I only: 



x = ± 




H = maxH(x), (6) 
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FIG. 1: (Color online). The shape of the base DC, H(x). 



where erf(z) = J* e~* dt. Solution (6) presents 

a self-affine one-parameter family of DCs parametrized 
by %. For a thin film, the DC describes the stationary 
profile of a drop with the amplitude % and zero contact 
angle. Owing to scaling (4), any DC can be expressed 
in terms of the base DC H{x) having H = 1, see Fig. 1. 
Thus, the profile of a DC and half its length Z obey: 



H(x) =HH (x/Vn 



1 



(7) 



The property of self-affinity is a necessary prerequi- 
site for the emergence of fractal structure described by 
Eq. (3), as we discuss below. As it follows from Eq. (7), 
DCs become narrower for smaller amplitudes - contrary 
to other examples of compactons where typically the 
width is amplitude independent [11]. In a more general 
situation, when for small h the functions f(h) and g(h) 
scale like f /g oc h 1 . no solutions are possible for 7 < — 2 
[13]. At 7 > 0, the solutions are soliton-like because their 
support is no more compact. The case 7 = results in 
constant Z. Compact solutions satisfying the requirement 
of Z — > as H — > exist in the range of — 2 < 7 < only. 
Thus, the considered above case 7 = — 1 is the only inte- 
ger index possessing self-affine compactons. 

To explore the stability of a DC, we introduce a small 
perturbation cx £(x)exp(At) of H(x), where A is the 
growth rate. By linearizing Eq. (3), we obtain 



= 0. 



(8) 



Assuming £(±Z) = 0, we multiply Eq. (8) by f" + H^ 1 ^ 
and integrate by parts to arrive at an integral relation 



a nonzero volume and can be destabilizing, if nonlinear 
corrections are retained. 

Thus, although a DC is stable with respect to per- 
turbations of zero volume, the instability is possible if 
the volume of the DC is changed, as confirmed by our 
numerical simulations of Eq. (3). We detect the tempo- 
ral decay even for finite-amplitude perturbations of zero 
volume. For the perturbations changing the volume but 
keeping constant the length of the DC, we find a breakup 
of the DC with the emergence of a complex structure. 

Evolutionary problem. — We now study the formation 
of a fractal, hierarchical structure of DCs, illustrated in 
Fig. 2. We discretize Eq. (3) in the computation do- 
main < x < d with periodic boundary conditions, with 
the number of nodes N = 1000 and apply the Newton- 
Kantorovich method [8]. We choose a distorted uniform 
profile h(x,t = 0) = 1 + acos(2nx/d) with a = 0.1 as an 
initial condition. Results of computations are presented 
in Figs. 2 and 3. There we also compare the numer- 
ically obtained profile h(x) having local maxima hm\ 
n = 1, 2, . . . with the DC profiles with U = h { m\ DC (n) 
below, what indicates that the initial state develops into 
a hierarchical structure of DCs of different amplitudes. 

This fact allows us to increase the efficiency of the nu- 
merics significantly: because after their formation the 
DCs remain constant in their bulk, we exclude these do- 
mains from numerical simulations and impose the corre- 
sponding boundary conditions for the still evolving do- 
mains between the formed DCs. Thus, while proceeding 
to smaller DCs, we can considerably refine the mesh and 
also increase the time step. Therefore, we not only re- 
solve high-order DCs with the accuracy consistent with 
that at previous stages, but also maintain the computa- 
tional efficiency. This strategy provides reliable results 
up to n = 4. 

The observed structure along with the property of self- 
affinity suggests that the formation of higher-order DCs 
never stops and the dry spots between DCs, form a fractal 
reminiscent of the Cantor set. Here, a DC plays a role 
of a primitive element, mediating the fractalization. To 
characterize properties of this fractal quantitatively, we 
plot in Fig. 4(a) the variation of L n , the distance between 
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dx, 



(9) 

which is closely related to the variational principle for 
Eq. (1) [14] and the fact that the Lyapunov functional has 
a local minimum on the DC. As H > 0, both the integrals 
in Eq. (9) are non-negative and the perturbations are 
nongrowing, A < 0. This result, however, does not guar- 
antee against the instability, as there exist two modes of 
neutral stability, A = 0, satisfying £(±Z) = 0. One mode, 
^J ' = H 1 , reflects translation invariance and cannot give 
rise to instability. Another mode, £2 = H — xff' /2, has 




FIG. 2: Fragment of the evolution of the field illustrating 
hierarchical formation of droplets. Notice logarithmic scales 
of the time and the field. 
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FIG. 3: (Color online). Evolution of the film profile for d = 
10. Panels (b) and (c) are zoomed in fragments of panel (a). 
Lines represent numerical results according to Eq. (3), circles 
show the profiles of corresponding DCs as in Eq. (6). 



the neighboring DCs of n-th and (n — l)-th orders, versus 
the base 2l n of DC'™- 1 . The numerical results for different 
d fit well a power law: 



L n fna(2l n f, aw 0.2, (3 ps 1.25. 



(10) 



Note that deviations from this law for the points related 
to the biggest DCs stem from the initial condition. On 
the other hand, for higher orders the self-similarity of the 
formation of DCs is evident from Fig. 4(a). 

Because j3 > 1 in Eq. (10), with the increase in n the 
ratio L n /l n diminishes implying that the smaller daugh- 
ter DCs tend to occupy the whole space between their 
bigger parent DCs. The fraction of dry spots tends 
to zero and, therefore, the fractal dimension of this 
set equals zero. Furthermore, for large n we can ne- 
glect the distance between DC (n) and DC (n+1) and put 
L n ps 2l n+ \. As a result, Eq. (10) entails a remarkable 
supercxponential scaling of l n with n: 



log(U cx/3 n log(Z ). 



(11) 




FIG. 4: (Color online), (a) The distance L n between two 
neighboring DC (n) and DC 1 "" 1 ' versus the base 2l n of DC (n) . 
Squares and circles are numerical results for d = 8 and d = 10. 
Dotted line is a fit, Eq. (10). (b) Mapping log Gfc+i(log Gk) 
calculated in the framework of Eq. (14) (circles); dotted line 
corresponds to the asymptotic law, Eq. (15). 



Self-similar solution. — To shed light on the hierarchi- 
cal formation of DCs and to alternatively support the 
conclusions about the fractal dimension and the super- 
exponential scaling, we construct self-similar solutions, 
which originate from the rescaling property, Eq. (4). By 
seeking the solution of Eq. (3) in the form 

h = t- 1 G{r)), r] = xVi, (12) 

we arrive at an ordinary differential equation for G[rf): 

77G'-2G + 2(G 2 G" + G 3 G'")' = 0, (13) 

where primes stand for d/drj. Numerical solutions of 
Eq. (13) with various initial conditions all demonstrate 
a qualitatively similar behavior of G(rj), which displays 
an infinite number of oscillations of increasing amplitude. 
Two numerical solutions corresponding to different initial 
conditions arc shown in Fig. 5. 

For large G, where we can estimate d/drj ~ e 1 / 2 with 
e ~ G _1 <C 1, the first two terms in Eq. (13) become 
negligible in comparison with the last two terms. In this 
limit, Eq. (13) is reduced to Eq. (5) with G{rf) instead 
of H(x). Therefore, G{rf) can be approximated by the 
solution for a DC [see Eq. (6) and the inset in Fig. 5] ev- 
erywhere except for its tails, where G is no longer large. 
As a result, G(rf) looks like a sequence of DCs with su- 
perexponentially growing amplitudes Gk ~ exp[^4 fc ] and 
widths Ai]k = 2v / 7rGfc; 
large k satisfy r]k ps y/irGk, 

To specify the superexponential growth of Gk with fc, 
we construct a mapping Gk —> Gk+i valid for large G. 
In the range of \rj — %| < \J irGk , G ps GkH(xk) with 
x,k = (77— ijk) /^JGk- To bridge the solution for DC^ 
with that for DC^ k+1 \ we substitute a representation 

G = e~ 2 C(y), V = (v-Vk- V^G k )e, s = Gfe 1/6 into 
Eq. (13) and neglect the terms oc e, what yields 

y oC' + 2(C 2 C' + C 3 C"')' = o. (14) 

Here, primes denote d/dy and yo = rjk/ WGk+V^ ^ 2-^7?. 
As Eq. (14) admits no analytical solution, we solved it 
numerically, with the initial condition at y = — £ _2 (v / 7r — 
x): d p C/dy p = e 2 ^-^d p H/dx p with p = 0,1,2,3 and 
d p H/dx p evaluated via Eq. (6) at H = 1, which ensures 
the matching with the decaying tail of DC^ . To perform 
the matching with the growing tail of DC (fc+1) at y > 1, 
we take into account that H" = — ln H — 1 [cf. Eq. (5)] 
and obtain 



the positions r\k of maxima for 
see markers in Fig. 5. 



G 



fe+i 



= Gi/ 3 Cexp(C" + l). 



By determining Q, we end up with the transformation of 
Gk — > Gk+i [see Fig. 4(b)], which fits well a power law 



G 



fc+i 



40 G 



2.83 



(15) 



Equation (15) shows a supercxponential growth for Gk 
with k, which is expected as required by the self-affinity 
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FIG. 5: (Color online). Function G(n) in solution (12). Initial 
conditions at small r\ are G w 0.5r? 2 (solid line) and G « 
1 — O.I77 2 (dashed line). Squares and triangles show the local 
maxima Gk of G(n); for large k, the maxima Gk approach the 
law G = rj 2 /it (dotted line). Inset: a comparison of a piece of 
G(rf) with a single DC (circles), Eq. (6). 

and the similar behavior for the lengths, see Eq. (11). As 
the amplitude oc I 2 [cf. Eq. (7)], the exponent 2.83 in 
Eq. (15) is in reasonable agreement with 2/3 in Eq. (10), 
obtained within the evolutionary problem. The fact that 
the correspondence is not perfect is not surprising as 
Eq. (15) is the asymptote of extremely large t (i.e. large 
k), while Eq. (10) is a fit obtained for the early stage of 
the evolution (small k). Nevertheless, we see that the 
self-similar solution is closely related to the hierarchical 
structure of DCs described by the evolutionary problem. 

Finally, we stress that the relation between the self- 
similar solution and the spatially periodic solution as in 
Fig. 3 is not simple. The whole structure of successive 
DC-like solutions, h(x,t), obtained via Eq. (13), moves 
with the time toward the point x = 0, whereas DCs 
H(x) 1 which are born as a result of evolution according 
to Eq. (3), are stationary objects. However, the long- 
time evolutions of both these solutions show the sim- 
ilar displacement of the gaps between DCs by higher- 
order DCs. This argument becomes transparent, if we 
observe the self-similar solution "stroboscopically" . Let 
us consider a self-similar solution at moments of time 
tk = Vk/ x o- The corresponding field profile (12) describes 
the formation of DCs up to the k-th order in the domain 
< x < xo+Arjk/y/tk with DC^- 1 centered at x = xq. As 
the growth of Gk with k is superexponential, the highest- 
order DC dominates the pattern, what ensures that the 
fractal made of the dry spots has zero dimension. 

Conclusions. — We have applied the concept of dissipa- 
tive compactons to the evolution of a thin film within 
a framework of the generalized one-dimensional Cahn- 
Hilliard equation. We have shown that as a result of a 
rupture, the thin film evolves into a hierarchical struc- 
ture of drops, which can be represented by dissipative 
compactons of different scales. By efficiently solving the 



amplitude equation and, alternatively, by constructing a 
self-similar solution, we show that this structure of DCs 
is a fractal, characterized by superexponentially decreas- 
ing amplitudes and lengths of smaller droplets, and thus 
having zero dimension. The dissipative compacton is a 
primitive element mediating the fractal structure com- 
prising the dry spots between the compactons. It should 
be also noted, that a number of effects, such as inter- 
molecular interaction between liquid and solid, contact 
angle dynamics, evaporation, etc. become of crucial im- 
portance, when the free fluid surface touches the solid. 
An extension of the theory above that includes these ef- 
fects remains a challenge. 
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